Astronomy & Astrophysics manuscript no. NaX'solarabund © ESO 2010 

May 25, 2010 



Highly Ionized Sodium X-ray Line Emission from the Solar Corona 

and the Abundance of Sodium 

K. J. H. Phillips 1 , K. M. Aggarwal 2 , E. Landi 3 , and F. P. Keenan 2 

1 UCL-Mullard Space Science Laboratory, Holmbury St Mary, Dorking, Surrey RH5 6NT, United Kingdom 
e-mail: kjhp@mssl.ucl .ac.uk 

2 Astrophysics Research Centre, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 INN, N. Ireland, 
United Kingdom 
e-mail: K.Aggarwal@qub.ac.uk, F.Keenan@qub.ac.uk 

3 Naval Research Laboratory, Washington, D.C. 20375-5320, U.S.A. 
e-mail: enrico . landi@nrl . navy . mil 

Received /Accepted 

ABSTRACT 

Context. The Nax X-ray lines between 10.9 and 1 1.2 A have attracted little attention but are of interest since they enable an estimate 
of the coronal abundance of Na to be made. This is of great interest in the continuing debate on the nature of the FIP (first ionization 
potential) effect. 

Aims. Observations of the Na x lines with the Solar Maximum Mission Flat Crystal Spectrometer and a rocket-borne X-ray spectrom- 
eter are used to measure the Na/Ne abundance ratio, i.e. the ratio of an element with very low FIP to one with high FIP. 
Methods. New atomic data are used to generate synthetic spectra which are compared with the observations, with temperature and 
the Na/Ne abundance ratio as free parameters. 

Results. Temperature estimates from the observations indicate that the line emission is principally from non-flaring active regions, 
and that the Na/Ne abundance ratio is 0.07 ± 50%. 

Conclusions. The Na/Ne abundance ratio is close to a coronal value for which the abundances of low-FIP elements (FIP < 10 eV) 
are enhanced by a factor of 3 to 4 over those found in the photosphere. For low-temperature (T e < 1.5 MK) spectra, the presence 
of Fe xvn lines requires that either a higher-temperature component is present or a revision of ionization or recombination rates is 
needed. 
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£T ! 1. Introduction 2 P 3/2 in the extreme ultraviolet, at 681.7 A and 694.3 A, and the 

i — ' „ ,. , „ ,, , , .... ,, He-like Na (Nax) line 1*2 .s 3 5 , i - ls2p 3 P? at 1111.76 A seen 

I ■ Sodium has lew prominent or unblended emission lines in the . „,.„,-. „ IT . „„ i„ ,, , A ~„^ n -„^7k j 

rfs , ,. • , . j v . •„ a u i ii m SOHO SUMER spectra (Curdt etal., 2000, 2001), and reso- 

<- ' ultraviolet and X-ray spectra emitted by solar coronal plasmas. .. ,„ ... ,., . 777 — Tr n ,x7 — T^ \, r^ 

■*-+■ ■ TT , . .. j... , , , .. , , nance lines ol H-like (Na xi) and He-like Na (Na x) in the solt X- 

NJ . , Hence determinations ol its abundance relative to other elements ,. . . . . XT ,. ' v 

. „ .. „ , ... . ray range. X-ray observations of the Nax and Naxi X-ray emis- 

l/^ rely on measurements ol emission line fluxes that are either weak ... , ,, . .... . , „ . . , 

^— ; . , , ., „ . . , , sion lines should in principle be easier than lor the equivalent 

Z^J or in crowded spectral regions. Yet its abundance is important . .. i i i ■ i j r xt • i • i 

0„ ,. . r „ , ,, . t-tt-. rc ii. ii potassium lines as the photosphenc abundance of Na is higher 

. lor discussions ol the so-called FIP effect, whereby coronal el- f „ , , , ^ r „ i . r . , ^ . i »„„„. , , 

. , , ,.„ c ,. , . J , . by a factor ol 16 than K (Asplundet al., 2009), but neverthe- 

^ ement abundances differ from corresponding photosphenc val- , . — r __ — rr - , . T _ ,. 

• • • , ,. , ., ., . ,. . . . . *. , /T , Tri s • less observations remain scarce. The Na xi Ly-a and Ly-/j lines 

^> ues depending on whether their first ionization potential (FIP) is » » • J -^ 

•^ ■ neater or less than about 10 eV JMeTel [1985; Fel dman et all are at 10 ' 02 A and 8 ' 46 A res P ectlvel y' and the Na x 1 , 2 - 1 S 2l 

X ■ 119929) . This is because the FIP of Na is only 5.14 eV, less il = Si P) lln f are near 11 A. The Ly-/3 line was observed by 

b ! than any other element common in the Sun apart from potas- I Walker et al. | dl974j) with crystal spectrometers on the OV1-17 

sium, which has an FIP of 4.34 eV. Determinations of the coro- spacecraft, arid the Nax lines with rocket-borne crystal spec- 

nal abundance of either Na or K can therefore probe the pos- trometers by |Parkinson| <tl975|). Only a single scan of the Flat 

sibility that elements with small FIP are enhanced by amounts Cl T stal Spectrometer (FCS) on Solar Maximum Mission (SMM) 

that depend on the magnit ude of the FIP, as i s sugge sted by the in its 9 "y ear lifetime was made of the Na x lines during a flare 

electric current model of iHenoux & S"oTr^vl fl99% to explain ln 198 ° (Mips et alj, U.982). Previous Na abundance analyses 

the FIP effect. Ob servational evidence for this has been pro- of the <5Vi-77 data gave log A(Na) = 6.26 (|Walker et alj.|1974h 

vided by Feldman ( 1993), and to some extent by observations and of the Parkinson rocket data log A(Na) = 6.73 (Parkinson, 

of He-like K X-ray hnes seen in numerous flares and non-flare 1975 ) ( on a scale lo § A ( H ) = 12 )- The Na x and Na XI X " r ay lines 

periods between 2002 and 2005 with the RESIK spectrometer have incidentally been observed in radiation from a laboratory 

(ISvlwesteretaQ l20T(ih . Coronal sodium emission lines include (Z-pmch) device by| Burkhalteret al. |(|1990|). 

the Li-like Na(Narx) resonance lines ls 2 2s 2 Si /2 -ls 2 2p 2 Pi /2 , Only rather approximate atomic data were available at the 

time of the Walker et al. ( 1974) and Parkinson ( 1975) analyses. 
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needed for the evaluation of the abundance of sodium from X- 
ray line fluxes. The chianti database and software package in the 
Interactive Data Language (I DL) Solar SoftWare system, since 
its inception in the late 1990s dDere et al.Lll997ll2009h . has also 
simplified analyses of X-ray and ultraviolet spectral data consid- 
erably. Excitation data for the Na x lines presently included in 
chianti are interpolated values from other elements, but more re- 
cently Aggarwal et al. (2009) have calculated collisional excita- 
tion data specifically for He-like sodium with the close-coupling 
/^-matrix code. Also, more refined ionization fractions assuming 
coronal ionization e quilibrium have recently been calculated by 
lBrvansetal.ld2009t) . 

The availability of new atomic data has inspired a fresh anal- 
ysis, reported here, of the Na x X-ray lines seen in the SMM 
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FCS spectral scan and the spectra obtained by Parkinso: 
X-ray lines of Na x were identified by Phillips et 



1975) 



1982), 



with m ore detailed analysis of nearby lines by Landi & Phillips 
(2005). A slight revision of spectral line identifications is given 
here and a determination of the Na/Ne abundance ratio for this 
flare plasma, based on a remarkable coincidence of the contribu- 
tion functions of the Na x w line and that of the Ne ix Is 2 - 1 sAp 
(w4 or He-y) line blended with the Na x w line. The Na/Ne abun- 
dance ratio is examined in the light of other recent abundance 
determinations and the nature of the FIP effect in solar coronal 
plasmas. 



2. Observations 



The SMM FCS observations dPhillips et all 119821) were made 
during the decay stage of a GOES Ml. 5 class flare on 1980 
August 25 with channel 2 of seven channels making up the in- 
strument. Scanning flat crystals were mounted on a rotatable 
shaft, and radiation incident on them was via a grid collima- 
tor giving a 14 arcsecond (FWHM) field-of-view. The high- 
precision Baldwin drive-encoder units used to rotate the drive 
led to very accurate wavelengths for observed spectral lines. The 
whole FCS could also be scanned spatially across an emitting 
region on the Sun. For the 1980 flare, a spatial scan was made 
over the flare emission near its onset to determine the location 
of the brightest point. The instrument was then pointed to this 
position, and a complete spectral scan taken. High-quality beryl 
(Be 3 Al 2 Si6 Ois, 2d - 15.96 A) was the diffracting crystal of 
channel 2, giving a spectral resolution of 0.0056 A at the wave- 
lengths of the Nax lines near 11 A. The scan took 17 minutes 
to accomplish, during which time there was a significant change 
in the emitting regi on. An ana l ysis o f the FCS data in this flare 
with chianti dLandi & Phillipsl r2005b took this into account, di- 
viding the scan into seven time bins. The Nax lines near 11 A 
lay at the end of time bin 1 and the start of time bin 2 of the 
scan, when the flare emission and temperature were still high. 
An emission measure analysis by Landi & Phillips (2005) us- 
ing only Fe lines (stages Fe xvn to Fe xxm) indicated a temper- 
ature of about 8 MK, which is reasonable for a Ml. 5 class flare 
a few minutes after its maximum emission. This analysis omit- 
ted emission lines of He-like ions such as Ne ix and Mg xi since 
they indicated a much lower temperature, of about 1-4 MK. The 
reason for this is that the emission measure appears to have a 
bimodal distribution, with a flare component (~ 8 MK) and a 
lower-temperature component attributable to the non-flaring host 
active region which gives rise to the emission from the He-like 
ions. The Nax emission lines appear to be in the latter cate- 
gory. Support for such emission measure distributions is pro- 
vided by analysis of Arxvn X-ray lines observed during flares 
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Fig. 1. X-ray spectrum in the 10.9-11.3 A range obtained with 
the SMM FCS during the decay of the 1980 August 25 flare 
(histogram) with theoretical spectra for temperatures 2 MK, 
3MK, 4 MK, and 5 MK (indicated in each panel) and an as- 
sumed Ne/Na abundance ratio equal to the coronal value of 0.07 
(Feldmanet al., 1992a). Line styles for the theoretical curves: 
continuous line = Na x lines alone (including weak Na ix satel- 
lites); dot-dash line = Nax and Neix lines; dotted line = with 
Fe xvn lines. The positions of the principal Na x lines w, y, and 
Z are indicated (the Na ix satellites j and k are within 0.004 A of 
line z), and the Ne ix w>4 line is at 1 1 .000 A. The other prominent 
lines are due to Fexvn, with observed wavelengths 11.023 A, 
11.133 A, and 11.253 A (see TableQ}. (A colour version is avail- 
able in the on-line journal. On-line journal version key: black 
histogram = observed spectrum; coloured curves are theoretical 
spectra, with code blue = Nax lines alone; green = Nax and 
Ne ix lines; red = with Fe xvn lines.) 



by the RESIK instrument dSvlwester et all 120081) : in that case, 
Arxvn line ratios were best fitted by a bimodal emission mea- 
sure distribution having temperatures of 4.5 MK and 16 MK. 
Also, analysis of broad- band data from the RHESSI instrument 
in its A0 attenuator state (Phillips et al., 2006) similarly indicates 
that there is a non-flaring component of emission, with temper- 
ature corresponding to the active region, as well as a hotter flare 
component. 

The FCS channel 2 spectral scan over the 10.9-1 1 .3 A range 
is shown in each of the four panels of FigureQ](histogram plot). 
In these plots, FCS count rates have been converted to absolute 
spectral irradiance units using pre-launch intensity calibration 
factors. The background emission is largely due to fluorescence 
of the crystal material with the solar continuum making a mi- 
nor contribution. The theoretical wavelengths of Na x lines are 
indicated. 

The rocket- borne crystal spectrometers described by 
iParkinsonl d 19751) were launched on a stabilised Skylark rocket 
on 197 1 November 30. The instruments viewed a non-flaring ac- 
tive region on the Sun through a multi-grid collimator with field 
of view equal to 3 arcminutes (FWHM). Two of the spectrome- 
ters scanned through the region of the Na x lines near 1 1 A with 
KAP (2d = 26.64 A) and gypsum (2d = 15.19 A) crystal. As 
the data are no longer in digital form, the spectra in the region of 
the Nax lines (10.9-1 1.4 A) from Figure 3 of lParkinsor] dl975h 
were hand-digitised for the purposes of this analysis. The photon 
count rate is higher in the KAP crystal scan and so the spectrum 
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Fig. 2. X-ray spectrum (histograms) in the 10.9-11.3 A range 
obtained with a rocket-borne KAP crystal spectrometer of a non- 
flaring active region in 1971 (Parkinson! 1 19751) . The spectrum 
is compared with synthetic spectra having temperatures 1 MK, 
1.5 MK, 2 MK, and 2.5 MK (indicated in each panel) and an 
assumed Ne/Na abundance ratio equal to a coronal value of 0.07. 
Na ix satellites are included in the synthetic spectra. Line styles 
are as for Figure Q] (A colour version is available in the on-line 
journal. On-line journal version key: black histogram = observed 
spectrum; coloured curves are theoretical spectra, with code blue 
= Nax lines alone; green = Nax and Neix lines; red = with 
Fexvn lines.) 



has high statistical quality. This scan is shown in Figure [2] The 
temperature is likely to be lower than the FCS scan in Figure Q] 
and the line feature at 1 1.152 A, if real, could be the Na ix satel- 
lites q and r; the theoretical position of line q is indicated. 

3. Atomic data and line identifications 

3.1. The spectrum in the 10.9-11.3 A region 

Na x lines in the X-ray region include a group of four lines with 
transitions of the type Is 2 - 1 s2l (I - s, p) a nd higher mem bers 
of the Is 2 - \snl series. Using the notation of lGabriell d 1 972h . the 
n — 2 lines are, in increasing order of wavelength, identified with 
letters w (transition Is 2 So - ls2p V P\, wavelength 1 1.003 A), x 
(ls 21 5 -ls2/? 3 P 1 ,11.080A),y(ls 21 5o-Fv2p 3 f , 2, H.083A), 
and z (Is 2 1 Sq - ls2s 3 S \, 11.192 A). The wavel engths, quoted 
from chianti, are from Ralchen ko et al.l d2008l) . Dielectronic 
satellites with transitions 1 s 2 n'l' - 1 s2pn'l' occur in great profu- 
sion in the neighbourhood of the w-z lines, but with only a few 
that are significant at lower temperatures, particularly satellites 
j and k (transitions 1 s 2 2p - 1 s2p 2 ) which blend with line z and 
satellites qandr {\s 2 2s-\s2s2p). Numerous satellites with very 
high n values occur near line w and converge on it; individually 
they are not important but their cumulative contribution is up to 
about 10% of line w. 

The Neix Is 2 1 Sq - ls4p l P\ (w4) line occurs very near the 
Nax w line, but how near was unclear from early atomic struc- 
ture calculations. Those of Vainshtein & Safronova ( 1978) gave 
the wavelength as 1 1 .025 A, longward of the Na x w line by 
an amount easily resolvable with t he FCS an d Parkin son instru- 
ments. Using the wavelengths of IWiese et alj (1196 9). Parkinson 
d 19751) was led to identify a line feature at 11.027 A as the 



Neix wA li n e. Ho wever, the more recent MCDF calculations of 
IChen et al . (2006) seem to have established that this line is at 
1 1 .000 A, just 0.003 A away from the Na x w line. In light of this, 
the feature at 1 1 .027 A is now identifiable with an Fe xvn line 
with transiti on 2s 2 2p 6 l S -2s2p 6 A p x Pi (wavelength 1 1 .023 A 
according to Landi & Phillips ( 2005)). More intense Fe xvn lines 
occur at 11.133 A and 11.253 A. Table[T]lists identifications and 
wavelengths of the principal lines in this region. 

Ne vm satellite lines with transitions 1 s 2 nl - IsnlAp are ex- 
pected to form weak features to the long-wavelength side of the 
Ne ix w4 line, but are outside the range of interest here: the most 
prominent are those in the 1 s 2 2p - \s2pAp array which are lo- 
cated at ~ 11.45 A. 



3.2. Contribution functions 

The contribution function G(T e ), defining the temperature range 
of significant emission of a line emitted by an ion (stage +m) of 
element X in a transitio n from an excited le vel j to the ground 
level, is defined by (e.g. lPhillips et alJ (12008b ') 



G(T e ) 



N(Xy i ) N(x +m ) N(X) N(H) 
N(X +m ) N(X) N(H)~ 



N e 



(1) 



where the number densities N are of the excited level of ion X| m , 
the ion X +m (all levels summed), the element X (all ionization 
stages m), and hydrogen (H), and N e is the electron density. The 
abundance of X relative to H is N(X)/N(U) - A(X), and for a 
coronal plasma N(H)/N e is 0.8. These functions can be deter- 
mined from chianti with the user supplying the line wavelength 
and a chosen set of ionization fractions and element abundances. 
Figure [3] shows the chianti G(T e ) functions for the lines of inter- 
est here, viz. Na x w, Ne rx w4, and the Fe xvn line at 1 1 . 129 A. 
The excitation data for the Na x lines used in these calculations 
are from interpolation of distorted-wave calculations for neigh- 
bouring ions, while t hose for the Neix lines are from the R- 
matrix calc ulations of Chen et al. (2006). The recent ionization 
fractions of iBrvans et al.l (120091) were u sed, as were the coro- 
nal abundances of Feld man et al.l (Il992al) . in which low-FIP el- 
ements have abundances enhanced over photospheric values by 
a factor of ~ 4 but high-FIP element abundances are equal to 
photospheric values (see §3.5). 

While the Fe xvn G(T e ) curve steeply rises at T e < 4 MK 
and falls at T e > 10 MK, reflecting the fractional abundance 
of the Fe +I6 ion, both the Nax and Neix lines have contribu- 
tion functions that are much more slowly decreasing with T e . 
There is a remarkable coincidence of the Na x and Ne ix curves 
for T e < 4 MK, and even at higher temperatures the Na x curve is 
consistently only a factor 3 above the Ne ix curve. This provides 
a means for determining the Na/Ne abundance if the two lines, 
which are practically unresolvable with the FCS and Parkinson 
spectrometers, are mostly emitted at temperatures < 4 MK. 



3.3. Theoretical Na X line fluxes 

The principal Nax X-ray lines in the 10.9-11.3 A range are 
mostly excited by electron collisions, and so collisional excita- 
tion rate coefficients are of considerable importance in synthe- 
sizing the spectra. As mentioned, the chianti data use interpo- 
lated distorted-wave calculations, but recent excitation data of 
Aggarwal et al. (2009) using the close-coupling /^-matrix code 
developed at Queen's University Belfast should be a significant 
improvement and were used here. The data are in the form of 
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Table 1. Principal lines in the 10.9-1 1.3 A region. 
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(a) Wavelengths from Landi & Phillips (2005). 
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Fig. 3. Contribution functions G(T e ) for lines in the 10.9-1 1 .2 A 
range: Naxw Is 2 1 S -ls2p 1 P l 1 1 .003 A; Ne ix w4 ls 2l S - 
\sAp A P y 11.000 A; Fexvu2s 2 2p 6 l S -2s 2 2p 5 5d l Pi 11.129 A. 
The contribution functions of the Na x and Ne ix lines are almost 
identical for T e < 3 MK. (A colour version is available in the 
on-line journal. On-line journal version key: blue curve = Na x 
w line; red = Ne ix w4 line; green = Fe xvn 1 1 . 129 A line.) 



temperature-averaged collision strengths Y(T e ) for the various 
transitions involved, and include auto-ionizing resonances which 
make a contribution to the Ts if near the excitation threshold. For 
the Naxw line, the excitation is mostly from the ground state 
Is 2 l So to the upper level ls2p P\, with spontaneous radiative 
transition from the upper level. However, for the x, y, and z lines 
in which the upper levels are triplets (l.$2.s 3 S\ or ls2p 3 f), 
other transitions are involved. Cascade transitions from higher 
levels 1 s nl where n > 3 are important, especially for lines x, y, 
and z, as are transitions from the upper levels of each line result- 
ing from recombination of the H-like stage of Na (Na +1 °). The 
new excitation data were entered into chianti data files so that 
Na x line fluxes could be obtained. While the differences in the 
x, y, and z line fluxes between the original interpolated data from 
chianti and the new /^-matrix results are rather modest (< 7 %), 
they are up to 22 % (T e = 3 MK) for the important w line. 

3.4. Na IX dielectronic satellites 

Dielectronic satellite line emission is an important contributor 
to X-ray spectra of He-like and H-like ions, particularly those 
with high atomic number Z. The chianti database currently has 
no dielectronic satellite excitation data for satellites due to Li- 
like Na which occur near the Na x X-ray lines, having transitions 



1 s 2 nl - 1 s2pnl (nl being the quantum numbers of the "spectator" 
electron in the transition: see Gabriel (1972)). Previous calcula- 
tions have used a variety of atomic codes to obtain these data; 
here we used an adaptation of the Cowa n pseudo-relativistic 
Hartree-Fock (HFR) code dCowanl 1 1981b for small personal 
computers designed by A. Kramida (2008, priv. comm.). 

The data concerned are the satellite (sat) wavelength and an 
intensity factor B sat . A satellite line is excited by dielectronic re- 
combination when an electron recombines with a He-like Na ion 
(Na +9 ) to produce a doubly excited state which then radiatively 
de-excites. The flux at the Earth (distance from Sun = 1 AU) of 
the satellite F^\ formed in this way is given by 



,+9 



F d 



iel 



2.07 x 10 



_ l6 N{Na +v )N e V A r A a exp{-E s Jk B T e ) 



4tt(AU) 2 A a + IA r 



n3/2 



(2) 



where A r and A a are the radiative and autoionization probabili- 
ties from the satellite's upper level respectively, E sai the energy 
of the upper level with respect to the ground level of the Na +9 
ion, V the flare emitting volume, and Icb is Boltzmann's constant. 
The term N(Na +9 )N e V can be rewritten 0.8/(Na +9 )A(Na)£M 
where the volume emission measure is EM — N 2 V and the frac- 
tion of Na +9 ions, JV(Na +9 )/2JV(Na +m ), is f(T e ). The intensity 
factor B sdt , defined by 



Ocat — 



A r A a 



(3) 



is calculated by the Cowan HFR code together with all the ra- 
diative and autoionization probabilities and the satellite wave- 
lengths. 

Input to the code is the satellite array and value of E sat . 
Single-electron radial functions are calculated by the code and 
combined with Slat er-Condon theory t o give energy levels. 
Previous experience dPhillips et al.Lll994l) has shown that choos- 
ing 100% for the scaling of the Slater parameters in the code 
gives good results for X-ray transitions, though there is gener- 
ally a small wavelength shift AA between the wavelength of the 
He-like ion line w and that of satellites with very high values of 
n for the spectator electron, which should converge on line w. In 
running the code for Naix satellites, we took 100% scaling of 
the Slater parameters and added AA = +0.001 A to the Cowan 
satellite wavelengths since this value of AA is needed to achieve 
agreement between the wavelengths of high-n satellites and the 
Na x w line. 

Table [2] lists data for a selection of satellites - wavelengths, 
Z? sat , and £ sat - from the calculations, which included satellites 
with spectator electrons having nl - 2s, 2p, 3s, 3p, 3d, etc., ex- 
tending up to nl = 6/j. We generated data for nl - lOp to deter- 
mine the wavelength shift AA. As was found by Gabriel d 19721) 
and others, satellites j and k (see Table [2} are the most intense, 
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Fig. 4. Maxwellian-averaged collision strengths T(T e ) as a func- 
tion of temperature for excitation from the Li-like Na (Na +8 ) ion 
ground level ls 2 2s 2 5 1/2 to the upper levels of satellites q, r, 
s, t, u, and v (transitions given in Table I3 calculated using the 
Flexible Atomic Code (FAC). (A colour version is available in 
the on-line journal, with different colours for the T curves for 
excitation to each upper level.) 



i.e. have the largest values of Z? sat . They are both within 0.004 A 
of the Nax z line at 11.192 A. This is in agreement with other 
calculations of satellite line wavelength s for elements having 
simila r atomic number, notably Mg (e.g. ISteenman-Clark et al.l 
(119841) '). For Ca (Z = 20), satellite k becomes resolved from line 
z, while for Fe (Z = 26), both j and k are resolved. 

Inner-shell excitation can also give rise to satellites in the 
ls 2 2s - \s2s2p array for solar flare densities (at much higher 
densities, inner-shell excitation gives rise to satellites in other 
arrays also). Again, the current version of chianti does not have 
atomic data relating to these satellites. As at least satellites q and 
r (see Tabled in the 1 s 2 2 s - 1 s2s2p array are expected to be im- 
portant, Maxwellian-averaged collision strengths Y(r<,) for the 
entire arra y were calculated using the Flexible Atomic Code of 
iGul (120031) . Assuming that the de-excitation is entirely the radia- 
tive transition to the ground level \s 2 2s 2 Si/2, the satellite line 
fluxes are then given by 

JV(Na +8 )JV e V 8.63 x !Q- 6 T(r e ) 
F i-s = , /tTT ,o x — exp (-E ex Jk B T e ) (4) 



4tt(AU) 2 



,1/2 



where Zs exc is the excitation energy of the satellite line. Figure|4] 
shows the calculated Y functions for excitation from the 
ls 2 2s 2 S 1/2 level to the levels indicated. 

3.5. Synthetic spectra 

A spectral synthesis program was written with the specific in- 
tention of matching the Na x and Ne ix lines, together with Na ix 
satellites, in the 10.9-1 1.3 A region as observed by the FCS and 
Parkinson spectrometer. It is assumed, as was found from the 
PsFSIK flare observations of Ar xvn line ratios dSvlwester et all 
2008) and analysis of broa d-band spec tra from RHESSI in its 
A0 attenuator state (Phillips et al., 2006J), that a two-component 
emission measure distribution describes the observed spectra, 
a low-temperature component (1-4 MK) appropriate for the 
Na x and Ne ix emission from the non-flaring active region and 
a higher-temperature component for the flare emission proper 



(FCS spectrum) or for a hotter part of the active region if present 
(Parkinson spectra): it is this component that is responsible for 
most of the Fe ion line emission in this region, including Fe xix 
and Fe xxm lines seen in the FCS spectrum an d which were fit- 
ted by an 8 MK component in the analysis by Landi & Phillips 
(120051) . The Fe xvn lines are emitted by plasma with T e > 3 MK 
(see Figure [3j> so contributions from both the active region and 
flare may be expected. The synthesis program computed spectra 
in a temperature grid from 0.8 MK to 5 MK. 

The Na x lines were in cluded in the synthesis program from 
fhe lAggarwal et al](l2009l) atomic data (§3.3) and the Na ix satel- 
lites, both those excited by dielectronic recombination and those 
formed by inner-shell excitation, from the data discussed in §3.4. 
The Ne ix w4 line was included using fluxes from chianti, as well 
as the much weaker intercombination line Is 2 So - ls4p 3 P\, 
unresolvable from the w4 line. Other lines were also included, 
most notably Fe xvn lines; lines of Fe xix and Fe xxm had neg- 
ligible fluxes in the temperature range considered. The chianti 
Fe xvn line wavelengths differ somewhat from the observed FCS 
wavelengths; as the FCS wavelengths are expected to be very 
precise (uncertainties < 2 mA in this regi on), th e FCS wave- 
lengths from an original analysis (Phillip s et all [1982) were 
used. Fluxes of the Na x w, x, y, and z lines were obtained from 
the collisional excitation data of Aggarwal et al.l (12009) which 
were inserted into files that could be read by the chianti software. 
The Na ix satellite fluxes were calculated from equations (2) and 
(4). Gaussian profiles were applied to spectral lines that were a 
convolution of the thermal Doppler broadening (FWHM width 
A/to) and the instrumental profile, defined by the crystal rocking 
curve (FWHM width equal to AA rc ), assumed to be Gaussian. 
The widths are given respectively by 



AA D = 1.665 



A /2^g7jo„ 

M ion 



1/2 



(5) 



where T mn is the temperature of the emitting ion, taken to be T e 
and M[ on the ion mass, and 



AA rr = 2d cos 8 Ad 



(6) 



where the crystal rocking curve is A6. A pre-launch measured 
value of A8 for the FCS observing the Nax lines of 100 arcsec- 
onds was taken for the FCS spectrum. The differing ion masses 
for the Na x, Ne ix, and Fe xvn lines were taken into account, 
though the FCS instrumental broadening dominates the con- 
volved line profile. For the Parkinson rocket spectra, the instru- 
ment profile was empirically taken to be 3 times the FCS rocking 
curve width to match the observed line profiles. As in the calcu- 
lation of contribution functio ns (§3.2), the N a and Ne ionization 
fractions were taken from Brv ans et al.l {2009). 

Element abundances w ere taken for an avera ge coronal 
plasma using values from iFeldman et al.l d!992al) (included 
as "coronal abundances" in chianti). An initial value of 
0.071 is thus taken for the sodium-to-neon abundance ratio - 
A(Na)/A(Ne) - but in our analysis it is a free parameter to be de- 
termined from line flux measurements. The Fe abundance like- 
wise is taken to be logA(Fe) = 8.10 (on a scale logA(H) = 12. 
This is supported by measurements of broad-band RHESSI flare 
spectra Philli ps et all d2006l) . indicating that of 27 flares, 19 
ha d estimated abundance s logA(Fe) that were within 20% of 
the lFeldman et alj d!992al) value. Altho ugh spatial var i ations i n 
coronal abundances have been noted (Feldma n et all Il992b). 
the chief variations are those during impulsive flares and above 
strong sunspot magnetic fields, when low-FIP elements appear 
to have photospheric abundances. 
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Table 2. Selected Naix satellite lines excited by dielectronic recombination: data from the Cowan HFR code. 



Wavelength (A) (fl) 


Line label (6) Transition 


Ssa, (S-') W 


E SM (Ry) 


11.043 


dl5 


s 2 3p 2 P l/2 -ls2p3pCP) 2 D 3/2 


3.33(13) 


73.5 


11.044 


dl3 


ls 2 3p 2 P 3/2 -U2p3p( [ P) 2 D 5/2 


5.77(13) 


73.5 


11.077 (rf) 


s 


ls 2 2s 2 S 1/2 - ls2p 2 ( 3 S) 2 P 3/2 


4.64(12) 


60.0 


11.078 w 


t 


ls 2 2s 2 S l/2 -ls2s2p( 3 S) 2 Pi/ 2 


2.90(12) 


60.0 


11.105 


n 


[s 2 2p 2 P m -\s2p 22 S m 


3.11(12) 


60.6 


11.108 


m 


s 2 2p 2 P 3j2 -\s2p 22 S m 


8.19(12) 


60.6 


11.152 (rf) 


1 


[s 2 2s 2 S i, 2 - \s2s2pCS) 2 P V2 


1.66(13) 


60.6 


11.155 w 


r 


ls 2 2s 2 S 1/2 -ls2s2p( i S) 2 P l/2 


9.61(12) 


60.6 


11.165 


b 


ls 2 2p 2 P 1/2 -ls2p 22 P 3/2 


3.28(11) 


60.6 


11.168 


d 


ls 2 2p 2 P 1/2 -ls2p 22 P ]/2 


7.12(10) 


60.6 


11.168 


a 


U 2 2p 2 P 3 , 2 -U2p 22 P 3l2 


2.13(12) 


60.6 


11.171 


c 


ls 2 2p 2 P 3/2 -U2p 22 P 1/2 


3.25(10) 


60.6 


11.188 


k 


s 2 2p 2 P l/2 -U2p 22 D m 


2.31(13) 


60.6 


11.191 


J 1 


[ S 2 2p 2 P 3/2 -U2p 22 D 5/2 


3.70(13) 


60.6 


11.192 


I 


s 2 2p 2 P 3/2 -U2p 22 D 3/2 


1.90(12) 


60.6 


11.286 


e 


s 2 2p 2 P V2 -\s2p 2A P, i2 


1.74(10) 


60.6 


11.298 <rf) 


u 


ls 2 2s 2 S 1/2 - ls2p 2 ( 3 S)*P V2 


3.45(9) 


60.0 


11.300 <d) 


V 


ls 2 2s 2 Si /2 -ls2s2p( 3 S) 4 P l/2 


6.64(8) 


60.0 


s from Cowan code + 


0.001 A; (b) Notation of Gabriel ( 1972); Belv-Dubau et al. ( 1979); (c) Numb 



powers of ten; (d) Lines also formed by inner-shell excitation. 



4. Results 

Figure [TJ shows the SMM FCS spectrum for the flare of 1980 
August 25. In this figure, the observed spectrum is compared 
with four synthetic spectra with T e = 2, 3, 4, 5 MK, all with 
the coronal value of A(Na)/A(Ne) = 0.071. In each case, the 
synthetic spectra were adjusted so that the blend of the Na x z 
line with the Naix j and k satellites (feature at 11.192 A) is fit- 
ted, leaving the blend of the Na ix w line and the Ne ix w4 line 
(feature at 11.003 A) free as well as the barely significant Nax 
y line at 11.083 A. At the temperatures shown, the Naix satel- 
lites j and k are very weak, the ratio of the sum of these two 
satellites to the Na x z line varying from 0.24 to 0.045 over the 
range T e — 2 MK to 5 MK. The satellites q and r are also weak, 
the feature at 1 1 . 15 A formed by them being unobserved by the 
FCS. In addition the ratio of blended Naix w to the Neix w4 
lines is only weakly T e -dependent (Figure [3). In fact, the only 
indication of temperature in the theoretical spectra in Figure [TJ 
is the presence of nearby Fexvn lines at 11.023 A, 11.129 A, 
and 1 1 .250 A. The best-fit temperature to the FCS spectrum is 
cl early T e — 3 MK. How ever, the emission measure analysis 
of Landi & Phillips (2005) shows that much of the emission of 
these and other more highly ionized Fe lines for this particu- 
lar stage of the August 25 flare occurs at a higher temperature, 
around T e ~ 8 MK, with He-like ion line emission such as Na x 
and Ne ix occurring at lower temperatures indicative of the host 
active region. The agreement of the FCS and theoretical spectra 
at T e — 3 MK for the Fe xvn lines in Figure [TJ may therefore be 
regarded as coincidental. 

As values of temperature over the range 2-5 MK only 
slightly affect the Na x and Ne ix line emission, the effect of the 
abundance ratio A(Na)/A(Ne) may be examined; this is shown in 
Figure [5] The FCS spectrum is compared with theoretical spec- 
tra, calculated at T e - 3 MK, for three values of A(Na)/A(Ne) : 
0.14, 0.07, 0.035, i.e. the coronal value o flFeldman etai1(ll992ah 
multiplied by 2, 1, and 0.5 respectively. Clearly the best match 
is for A(Na)/A(Ne) = 0.07. Ruling out the values of 0.14 and 
0.035 for this ratio (top and bottom panels of Figure [5) leads 
to an estimated precision that is approximately 50%. This esti- 
mate of the abundance ratio is significantly larger than the pho- 



tospheric abundance ratio, A(Na)/A(Ne) — 0.02 ( Asplund et al., 
2009), and suggests an enhancement of about 3 to 4 in the abun- 
dance of Na in the corona if the active region plasma is coronal 
in origin. 

For the Parkinson rocket spectra, the emission arises from 
a lower-temperature plasma indicative of a non-flaring active 
region. Figure [2] shows the KAP crystal scans compared with 
four theoretical spectra with T e = 1, 1.5, 2, and 2.5 MK and 
A(Na)/A(Ne) = 0.07. At T e < 2 MK, the ratio of the^ Nax w- 
Ne ix vv4 blend at 1 1 .002 A to the Na x z line at 1 1 . 1 9 1 A is sensi- 
tive to T e , largely through the contribution of the Na ix satellites 
j and k to Na x line z. The ratio of the sum of these satellites to 
Nax line z over the temperature range 1 MK to 2.5 MK varies 
from 2.5 to 0.14. Of the four theoretical spectra shown, the one 
with T e — 1 MK best fits the KAP spectrum. Although only just 
significant (total photon counts per 0.025 s equal to about 16), 
the feature at 1 1.15 A can be identified with the blend of Naix 
satellites q and r; the appearance of these satellites confirms the 
temperature of ~ 1 MK. This conclusion is also supported by 
the gypsum crystal scan though this h as low er statistical quality. 
In the original analysis of Parkinson (1975), an emission mea- 
sure distribution peaking at 3 MK was found, rather more than 
is found here with improved atomic data. The comparatively low 
temperature for this active region that we find is in keeping with 
the fact that solar activity was at a rather low level when the 
observations were taken. 

If the temperature of the active region is only ~ 1 MK, the 
Fexvn lines at 11.023 A, 11.129 A, and 11.250 A should no t 
be so strong. The ionization fractions of Brva ns et al.l {2009) 
indicate that there is a negligible fraction of Fe +16 ions for 
T e < 1.5 MK (Figure |3}. As with the FCS spectra, a higher- 
temperature component might be present. However, in that case 
the Na x and Ne x line emission would be much higher, the fit- 
ted temperature would be correspondingly higher, and disagree- 
ments in the fit to the Nax w-Neix w4 blend and the weak 
11.15 A line feature, formed by satellites q and r, would re- 
sult. There is a possibility, then, that the Fe ion fractions 
at T e < 1.5 MK are in error, and that there is in fact a non- 
negligible fraction of Fe +I6 ions at very low temperatures. 



K. J. H. Phillips et al.: Highly Ionized Sodium X-ray Line Emission from the Solar Corona and the Abundance of Sodium 




1 1.1 

Wavelength (A) 

Fig. 5. The FCS spectrum in the 10.9-11.3 A range during the 
flare of 1980 August 25 compared with theoretical spectra for 
T e — 3 MK and three values of the abundance ratio A(Na)/A(Ne) 
(indicated in each plot). Line styles for the theoretical curves as 
for FigureQ] Two of the Fe xvn lines are off-scale to show better 
the agreement of the theoretical curves with the observed Na x 
and Neix line emission. (A colour version is available in the 
on-line journal. On-line journal version key is the same as for 
FigureQ]) 



5. Conclusions 

In this work, recent collisional excitation results for Nax X- 
ray line emission have been applied to solar spectra for the first 
time, together with new atomic data for Na ix satellites, formed 
by both dielectronic recombination and inner-shell excitation. 
From these and other data from the chianti atomic code, syn- 
thetic spectra were calculated as a function of electron temper- 
ature T e . To date, the only observed spectra taken with high- 
resolution spectrometers in the 10.9-11.3 A range are from the 
SMM Flat Crystal Spectrometer (a single scan during the decay 
of a Ml .5 flare) and two rocket-borne scanning crystal spectrom- 
eters in 1971 viewing a non-flaring active region. Comparison of 
the FCS and calculated spectra in the 10.9-1 1.3 A range shows 
that the Nax X-ray lines with nearby Naix dielectronic satel- 
lites yield a Na/Ne abundance ratio through the blend of the 
Nax w (resonance) line with the Neix h4 (1 s 2 'So — ls4p l P\) 
line, these lines having almost identical contribution functions 
G(T e ) for the low temperatures considered here. The value ob- 
tained, A(Na)/A(Ne) = 0.07 + 50%, is the coronal value of 
Feld man et al] (Il992al) to within uncertainties. It is significantly 
higher than the photospheric value, 0.020, and suggests that 
sodium (like other low-FIP elements) is enhanced by a factor 
of 3 to 4 in the corona. The lower-temperature observations ob- 
tained by Parkinson ( 1975) are best fitted with theoretical spec- 
tra having T e = 1 MK, and are again consistent with a coronal 
Na/Ne abundance ratio. The presence of Fe xvn lines in these 
spectra is une xpected, the most re cent ionization equilibrium 
calculations of iBrvans et al] d2009l) indicating either that there 
is a higher-temperature component present or (more likely) that 



and advice on his spectra, as is Dr M. F. Gu for advice on the Flexible Atomic 
Code used in this work. 
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